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ABSTRACT 

We present deep Spitzer/IRS spectra for complete samples of 46 2Jy radio galaxies 
(0.05<z<0.7) and 19 3CRR FRII radio galaxies (z<0.1), and use the detection of polycyclic 
aromatic hydrocarbon (PAH) features to examine the incidence of contemporaneous star forma- 
tion and radio-loud AGN activity. Our analysis reveals PAH features in only a minority (30%) of 
the objects with good IRS spectra. Using the wealth of complementary data available for the 2Jy 
and 3CRR samples we make detailed comparisons between a range of star formation diagnostics: 
optical continuum spectroscopy, mid- to far-IR (MFIR) color, far-IR excess and PAH detection. 
There is good agreement between the various diagnostic techniques: most candidates identified 
to have star formation activity on the basis of PAH detection are also identified using at least 
two of the other techniques. We find that only 35% of the combined 2Jy and 3CRR sample show 
evidence for recent star formation activity (RSFA) at optical and/or MFIR wavelengths. This 
result argues strongly against the idea of a close link between starburst and powerful radio-loud 
AGN activity, reinforcing the view that, although a large fraction of powerful radio galaxies may 
be triggered in galaxy interactions, only a minority are triggered at the peaks of star formation 
activity in major, gas-rich mergers. However, we find that compact radio sources (D < 15 kpc) 
show a significantly higher incidence of RSFA (>75%) than their more extended counterparts 
(wl5 - 25%). We discuss this result in the context of a possible bias towards the selection of 
compact radio sources triggered in gas-rich environments. 

Subject headings: galaxies:active - infrared:galaxies 
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1. Introduction 

Extragalactic radio sources form an important 
subset of the Active Galactic Nuclei (AGN) popu- 
lation because their relativistic jets heat the Intcr- 
SteUar/Galactic Medium (ISM/IGM) and drive 
outflows on scale s ranging from th e kpc-scale nar- 
row line regions ( Holt et allbOOSl) . to the 100 kpc 
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1 Mpc hot X-ray haloes of the host galaxie s 



and galaxy clusters ()McNamara fc NulsenI 120071 ) 



Such feedback has been used to e xplain the coolin g 



problem for the hot X-ray gas (jBest et al.ll2005[) . 
as well as the hi gh end shape of the galaxy lumi- 
nosity function (jBenson et al.ll2003i) . and the cor- 
relations b etween black hole mass and host galax y 
properties ( Bower et al. 2006 : Croton et al. 2006[ ). 
Therefore, given that a large percentage of giant 
elliptic al galaxies may go through a radio-loud 
phase ( Best et all 2006[ ). the triggering of radio 
galaxies is a key issue for the understanding of the 
evolution of massive galaxies in general. However, 
there remain considerable uncertainties about how 
and when radio-loud AGN activity is triggered. 

Morphological and spectroscopic studies pro- 
vide strong evidence that powerful radio galaxies 
are triggered in galaxy mergers, although not nec- 
essarily at a single phase of a particular type of 
merger ("Hcckman et al. ll986tJ iRamos Almeida et al.l 



r2011b, a: .Tadhunter et al.ll201l[) . Given that hv 



drodynamic simulations of major, gas-rich merg- 
ers show that, as well as triggering AGN, such 
interactions are capa ble of triggering power- 
ful starburst activity ( Mihos fc Hernquisti 1996t 



Pi Matteo et al.l2007l:ICox et al.lbOOSH Johansso n et al 



20091 ) ■ the study of the level of recent star forma- 
tion activity (hereafter RSFAj^ in radio galaxies 
provides key information about the nature of the 
triggering events, as well as the timing of the AGN 
relative to any merger-induced starburst. 

The fact that radio-loud AGN are overwhelm- 
ingly hosted by early-type galaxies allows par- 
ticularly clean searches to be made for the signs 
of RSFA. Perhaps surprisingly, given the evidence 
that many of them have been involved in galaxy in- 
teracti ons ( Heckman et al.flQSdlRamos Almeida et al 
2011bl la[). only a minority of radio galaxies are 



found to show spect roscopic evidence for RSFA a t 
optical wavelengths (|Tadhunter et al.ll2002l lioul ) . 

However, it is possible that a substantial percent- 
age of any star formation activity is obscured at 
short wavelengths by circum-nuclear dust. More- 



^We define "recent star formation activity" to include all 
star formation activity that has occurred within 2 Gyr 
of the observation epoch, encompassing contemporaneous 
starbursts, continuous star formation and post-starburst 
stellar populations (e.g. post-starburst activity as detected 
in spectral synthesis modelling of optical spectra: Tad- 
hunter et al. 2005, Holt et al. 2007, Tadhunter et al. 2011). 



over, in some cases the optical continuum features 
that are characteristic of RSFA may be masked 
by AGN-related continuum emission components 
such as direct o r scattered AGN light, and neb- 
ular continuum (jTadhunter et al.l 120021) . Clearly, 
it is important to explore star formation diagnos- 
tics that are less sensitive to dust extinction and 
AGN-related continuum emission. Therefore we 
are undertaking a programme to investigate the 
degree of RSFA in radio galaxies using mid- to 
far-IR (MFIR) diagnostics. 

As a first step we used the Multi band Imag- 
ing P hotometer for Spitzer (MIPS: iRieke etal 



20041 ) to make deep MFIR photometry measure- 
ments of the complete 2Jy sample of 46 radio- 
loud AGN ( 0.05 <z < 0.7) ( Progr a m 20233: PI 
Tadhunter; iTadhunter et al.l l2007l : IPicken et al 
JOOS, gdoi hereafter T07, DOS, D09) as well as 
a sample of 19 nearby 3CRR radio-loud AGN 
from the Spitzer archive ( Dicken et al.ll2010i here- 



after DIO). The results show that [OIII] opti- 
cal emission line luminosity (ijoiiij) is strongly 
correlated with both the mid- (24/im) and far- 
IR (70/im) luminosities (L24/im and Lyo^m re- 
spectively). Since the [OIIIJA5007 emission from 
the NLR provides a good indication of the in- 
trinsic power of the illuminating AGN (e.g. 



Rawlings & Sau 


ndersll991; Tadhunter et al.ll998l 


Simpson! Il998l: 


LaMassa et al.ll2010: and discus- 



sion in D09), the correlations between MFIR lumi- 
nosity and [OIII] optical emission line luminosity 
provide strong empirical evidence to support AGN 
illumination as the dominant heating mechanism 
of the thermal MFIR emitting dust. We also used 
energetic arguments to demonstrate that, while 
much of the mid-IR emission is likely to be ra- 
diated by the warm dust in the torus, the NLR 
clouds are a plausible location for the cooler, far- 
IR emitting dust (D09). 

In addition, we found evidence for enhanced 
far-IR emission in the minority of radio galaxies 
that show evidence for RSFA activity at optical 
wavelengths. Our interpretation was that, while 
AGN illumination is the primary heating mecha- 
nism for both the warm (mid-IR emitting, 24/im) 
and cool (far-IR emitting, 70/im) dust in most 
powerful radio-loud AGN, heating by starbursts 
acts to substantially boost the 70/fm luminosity 
in the 20-30% of objects in the 2Jy sample with 
optical evidence for RSFA. 
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Unfortunately, apart from the most extreme ob- 
jects with the highest degrees of star formation ac- 
tivity, the far-IR excess does not necessarily pro- 
vide an accurate indication of the level of star for- 
mation activity in individual radio galaxies. This 
is because the degree of contamination of the far- 
IR continuum by emission from the NLR and 
circum-nuclear torus may vary substantially from 
object to object. An alternative is to use mid-IR 
spectroscopy to detect the strong Polycyclic Aro- 
matic Hydrocarbon (PAH) emission features at 
6.2, 7.7, 8.6 and 11.3/im. These well-defined emis- 
sion bands are considere d to provide an unambigu- 
ous s ignature of RSFA ( Diamond-Stanic fc Riekd 
20101) . and are much less susceptible to dust ex- 
tinction effects than optical continuum studies. 

Therefore we have undertaken a campaign 
to observe the complete 2Jy sample of 46 ob- 
jects (0.05 < z < 0.7), previously observed with 
Spitzer/MI PS, using th e Spitzer Infrared Spectro- 
graph flRSi lHouckeTal . 2004; program 50588: PI 
Tadhunter ) . In a similar manner to the work pre- 
sented in DIO, we also examine IRS spectra for a 
complete sample of 19 3CRR FRII radio galaxies 
that have, on average, lower redshifts and radio 
powers than the 2Jy sample. 

We aim to investigate the degree to which the 
presence and strength of the PAH emission fea- 
tures correlates with other diagnostics for identi- 
fying RSFA. In particular: is the detection of PAH 
features, and therefore RSFA by association, con- 
fined to the objects in which we already find evi- 
dence for such activity at optical wavelengths? Or 
is there a substantial population of radio galaxies 
in which the star formation activity is hidden by 
dust and/or masked by AGN- related continuum 
components at optical wavelengths? 

In this, the first paper related to our IRS pro- 
gram of spectral observations of powerful radio 
galaxies, we present an atlas of Spitzer IRS mid- 
IR spectra for the 2Jy and 3CRR samples, and 
discuss the results from the PAH emission analy- 
sis. In a second paper we will discuss the silicate 
features and fine structure lines detected in the 
spectra in the context of the unified schemes for 
AGN. We assume a cosmology with Ho = 71km 
s~iMpc-\rj„ = 0.27 and fix = 0.73 throughout 
this work. 



2. Sample Selection 

The primary sample consists of all 46 pow- 
erful radio galaxies and steep-spectrum quasars 
(F^ cx v- "',aii > 0.5) selected from the 2Jy 



sample of IWall fc Peacockl (|1985I ) with redshifts 
0.05 < z < 0.7, flux densities S'2.7ghz > 2 Jy 
and declinations S < 10°; the sample is complete 
based on these criterion. This sample is identi- 



cal to that presented in iTadhunter et al. (Il993h 



and Morganti et al.l (jl993l) except that the redshift 
and steep spectrum selection criteria have been 
applied, and the object PKS 0347-f05, which has 
since proved to fulfill the sam e selection criteria 
( di Serego-Alighieri et al.l [1994 ). has been added. 
The spectral index cut has been set to ensure that 
all the sources in the sample are dominated by 
steep spectrum lobe emission, while the lower red- 
shift limit has been set to ensure that these galax- 
ies are genuinely powerful sources. In the following 
we will refer to this sample as the 2Jy sample. 

The mid-IR spectra presented here complement 
a wealth of data that has been obtained for the 
2Jy sample over the last two decades. These 
i nclude: deep optical spectro s copic observations 



( Tadhunter et al. 1995 
2OO2I: iHolt et all I2OO7 



1998, 2002; Wills et al 



extensive observations 



at ra dio wavelengths (jMorganti et al.lll993lll997 . 
19991 : D0 8); complete deep optical ima ging from 
Gemini (jRamos Almeida et al. 2011bl la): deep 



Spitzer/MIPS mid- to far-IR photometric obser- 
vations (D08; detection rates 100% at 24/im and 
90% at 70/im). In addition, 85% of the sample 
have recently been observed with Chandra and/or 
XMM at X-ray wavelengths, including all objects 
with redshifts z < 0.2; and 78% of the sample have 
deep 2.2 micron ( K-band) near infrared imaging 
(jinskip et al.ll2010l) . The full sample of 46 objects 
includes a mixture of broad- line radio galaxies and 
radio-loud quasars (BLRG/Q: 35%), narrow-line 
radio galaxies (NLRG: 43%), and weak-line ra- 
dio galaxiefH (WLRG: 22%). In terms of radio 
morphological classification, the sample comprises 
72% FRII sources, 13% FRI sources, and 15% 
compact steep spectrum (CSS) /gigahertz peak 
spectrum (GPS) objects. The sample is presented 
in Table 



^ WLRGs are defined as having EW{[OIII]) <10 A, 
l lTadliunter et al.|[l99l) . 

^Additional information on tlicse objects and the 
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We also present IRS spectra for a complete 
sub-sample of 19 3C RR radio-loud AGN selected 
from the sample of Laing et al. ( 1983[ ) (see Ta- 
ble [5]). We have limited this sample to 3CRR 
objects with FRII radio morphologies and red- 
shifts z<0.1, leading to a sample which is complete 
in both Spitzer/MIPS detections (100% at 24/im 
and 89% at 70/im) and [OIII] A5007 emission line 
flux measurements (100%). The full sample of 
19 objects also includes a mixture of BLRG/Q 
(16%), NLRG (58%), WLRG (26%). Because the 
3CRR objects have lower radio powers and red- 
shifts on average than most of the 2Jy sample, 
they help to fill in the lower luminosity ends of the 
MFIR vs [OIII] correlations (DIG). In the follow- 
ing discussion we will refer to this sample as the 
3CRR sample. Note that, although two objects in 
the sample (3C277.3, 3C293) have uncertain ra- 
dio morphological classifications, and cannot be 
confidently characterized as either FRI or FRII 
types, they are included here for completeness. 
Also, the compact steep spectrum radio source 
3C305 is included because it presents a miniature 
FRII radio morphology, with distinct hotspots at 
the ends of its jets, even if its radio structure 
as a whole is highly distorted because of strong 
jet-cloud interactions (jHeckman et al.lll982l) . The 
Spitzer/MIPS flux data and associated errors are 
presented in DIO. The [OIII] fluxes were obtained 
from published deep optical spectra at both high 
and l ow resolution taken us ing Dolores on the 
TNG (|Buttiglione et al.ll2009t ). except for DA240, 
4C73.08, 3C321 and 3C445 (see DIO for details). 
Note that two objects overlap between the 3CRR 
and 2Jy samples (3C403, 3C445). 

3. Observations and data reduction 

The mid-IR spectra of 35 out of the 46 ob- 
jects in the 2Jy sample were obtained in a ded- 
icated campaign of Spitzer IRS observations (pro- 
gram 50558; RI. Tadhunter) between July 2008 
and March 2009. These data were taken in the 
short low (SL) and long low (LL) resolution star- 
ing modes, covering 5.2 to 14.5/im and 14.0 to 
38.0/im respectively. The resolution for the ob- 
servations was Rsa60-127 and the slit widths were 
3.6" and 10.5" for the SL and LL observations 



more extended 2Jy sa mple 
|http: / /2jy.extragalactic.info| 



can be found at 



respectively. The spectra of 8 additional objects 
in the 2Jy sample, taken between June 2004 and 
May 2006, were obtained from the Spitzer archive. 
These data came from various campaigns under 
several different Pis and, consequently, have vary- 
ing integration times and observing modes. De- 
tails of the observations can be found in Table [T] 

Out of the 46 objects in the 2Jy sample, 43 
(93%) were successfully detected by Spitzer/IRS. 
PKS2135-14 was not observed because the Spitzer 
observatory ran out of cryogenic coolant before the 
observations could be made, while no observations 
of PKS2211-17 were taken because its mid-IR flux 
level (0.5 mJy at 24/im) was deemed too low to de- 
tect in reasonable integrations with Spitzer IRS. 
Observations of a further object (PKS1648-1-05) 
were available in the Spitzer archive, however, 
the integration times were an order of magnitude 
lower than those used for objects of similar mid-IR 
flux in our program. Therefore, this source was 
not detected with the IRS instrument. The SL 
spectra of an additional three objects (PKS0034- 
01, PKS0035-02, PKS1602-I-01) obtained from the 
Spitzer archive potentially suffer from higher flux 
calibration uncertainties, due to the effect of sat- 
uration of part of the detector used for the short 
wavelength range observations. This problem oc- 
curs because the IRS detectors encompass not 
only the spectral 2D images, but also the aper- 
tures used for the target acquisition imaging. If 
the ambient infrared background is high (above 
25 mjy/sr) then, for long integration times, the 
acquisition images are likely to saturate, causing 
the whole detector to perform irregularly. 

The majority of observations (40/43) were 
taken in Staring Mode, in which the target is ob- 
served in several cycles of set integration time and, 
for each cycle, the IRS instrument nods the object 
between two positions on the detector. For every 
cycle the nod positions were subtracted (A from 
B, and B from A) in order to perform background 
subtraction on the image, as well as to counteract 
intrinsic sensitivity variations across the detec- 
tor. The two resulting spectra were combined to 
produce the final spectrum. 

Spectra for three 2Jy sources, obtained from 
the Spitzer archive, were observed in Mapping 
Mode in which the spacecraft makes a pre- 
determined number of pointings of set integra- 
tion time mapped over a target. The pattern can 
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Table 1 

The 2Jy Sample - Spitzer Observational Data 



Name 


Other name 


z 


Mode 


Obs. date 




cycles 


tint {^bL2J 


cycles 


+ /T T 1 \ 


cycles 


*int (LL2) 


cycles 


XbLi 


0023-26 




0.322 


S 


,Ian-09 


60.0 


8 


60 


8 


120 


4 


120 


4 




0034-01 


3C15 


0.073 


S 


Jan-06 


240 


2 


240 


3 


120 


3 


120 


2 


0.8 


0035-02 


3C17 


0.220 


s 


Jaii-06 


240 


2 


240 


3 


120 


3 


120 


2 




0038+09 


3C18 


0.188 


s 


Jan-09 


14 


16 


14 


16 


30 


4 


30 


4 




0039-44 




0.346 


s 


Dec-08 


60 


4 


60 


4 


30 


4 


30 


4 


1.1 


0043-42 




0.116 


s 


Doc-08 


60 


4 


60 


4 


120 


2 


120 


2 




0105-16 


3C32 


0.400 


s 


,Tan-09 


60 


4 


60 


4 


120 


2 


120 


2 




0117-15 


3C38 


0.565 


s 


,Tan-09 


60 


8 


60 


8 


120 


4 


120 


4 




0213-13 


3C62 


0.147 


s 


Jan-09 


60 


2 


60 


2 


30 


4 


30 


4 


1.1 


0235-19 


OD-159 


0.620 


s 


Oct-08 


60 


4 


60 


4 


120 


2 


120 


2 




0252-71 




0.566 


s 


Oct-08 


240 


2 


240 


2 


120 


4 


120 


4 




0347+05 




0.339 


s 


Oct-08 


60 


8 


60 


8 


120 


4 


120 


4 




0349-27 




0.066 


s 


Oct-08 


60 


4 


60 


4 


120 


2 


120 


2 




0404+03 


3C105 


0.089 


M 


Sep-05 


14 


1 


14 


1 


14 


1 


14 


1 




0409-75 




0.693 


s 


Dec-08 


240 


2 


240 


2 


120 


4 


120 


4 




0442-28 




0.147 


s 


Nov-08 


60 


4 


60 


4 


30 


4 


30 


4 




0620-52 




0.051 


s 


Dcc-08 


240 


2 


240 


2 


120 


3 


120 


3 


1.2 


0625-35 


OH-342 


0.055 


s 


Dec-08 


60 


2 


60 


2 


30 


4 


30 


4 


1.1 


0625-53 




0.054 


s 


Dec-08 


240 


2 


240 


2 


120 


4 


120 


4 




0806-10 


3C195 


0.110 


s 


Dec-08 


14 


4 


6 


2 


6 


2 


6 


2 


1.1 


0859-25 




0.305 


s 


Jan-09 


60 


4 


60 


4 


120 


2 


120 


2 


1.3 


0915-11 


Hydra A 


0.054 


s 


Dec-05 


60 


9 


60 


11 


30 


4 


30 


4 


3.0 


0945+07 


3C227 


0.086 


M 


May-06 


14 


1 


14 


1 


14 


1 


14 


1 


1.2 


1136-13 




0.554 


s 


Jul-08 


60 


4 


60 


4 


120 


2 


120 


2 




1151-34 




0.258 


s 


Fob-09 


60 


2 


60 


2 


30 


4 


30 


4 




1306-09 




0.464 


s 


Aug-08 


14 


32 


14 


32 


120 


3 


120 


3 




1355-41 




0.313 


s 


Mar-09 


60 


2 


60 


2 


30 


4 


30 


4 


1.1 


1547-79 




0.483 


s 


Oct-08 


60 


4 


60 


4 


120 


2 


120 


2 


1.1 


1559+02 


3C327 


0.104 


s 


Oct-08 


14 


4 


6 


2 


6 


2 


6 


2 


1.1 


1602+01 


3C327.1 


0.462 


s 


Mar-06 


240 


2 


240 


3 


120 


3 


120 


2 




1648+05 


Here A 


0.154 


s 






















1733-56 




0.098 


s 


Oct-08 


60 


2 


60 


2 


30 


2 


30 


2 




1814-63 




0.063 


s 


Nov-08 


60 


2 


60 


2 


30 


2 


20 


2 


1.1 


1839-48 




0.112 


s 


Nov-08 


60 


6 


60 


6 


120 


4 


120 


4 




1932-46 




0.231 


s 


Nov-08 


60 


8 


60 


8 


120 


4 


120 


4 


0.7 


1934-63 




0.183 


s 


Nov-08 


60 


4 


60 


4 


30 


4 


30 


4 




1938-15 




0.452 


s 


Nov-08 


14 


16 


14 


16 


120 


2 


120 


2 


1.8 


1949+02 


3C403 


0.059 


M 


Oct-05 


14 


1 


14 


1 


14 


1 


14 


1 




1954-55 




0.060 


s 


Nov-08 


60 


6 


60 


6 


120 


4 


120 


4 




2135-14 




0.200 


s 






















2135-20 


OX-258 


0.635 


s 


Nov-08 


14 


28 


14 


28 


120 


4 


120 


4 




2211-17 


3C444 


0.153 


s 






















2221-02 


3C445 


0.057 


s 


Jun-04 


30 


4 


60 


2 


14 


2 


14 


2 


1.2 


2250-41 




0.310 


s 


Nov-08 


60 


4 


60 


4 


120 


3 


120 


3 




2314+03 


3C459 


0.220 


s 


Jan-09 


14 


16 


14 


16 


30 


2 


30 


2 




2356-61 




0.096 


s 


Nov-08 


60 


4 


60 


4 


30 


4 


30 


4 





Note. — Table presenting the basic parameters for the 2Jy sample as well as details of the Spitzer/IRS observations. Column 
4 presents the IRS observing mode: S — Staring mode; M - Mapping mode. Integration time in seconds and cycles for the IRS 
observations are presented in columns 6-13. SLl: 5.2— 8.7/i.m; SL2: 7.4— 14.5/i.m; LLl: 14— 21.3/i.m; LL2:19.5— 38. Column 14 presents 
the multiplication factor applied to the SL spectrum in order to match the LL spectrum where applicable. 
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Table 2 

The 3CRR Sample - Spitzer Observational Data 



Name 


z 


Mode 


Obs. date ti 


nt(SLl) cycles ti 


„t(SL2) 


cycles ti 


„t(LLl) 


cycles ti 


^t(LL2) 


cycles 


xSL 


3C33 


0.060 


S 


J an- 05 


60 


2 


60 


2 


120 


2 


120 


1 


1.0 


3C35 


0.067 


M 


Jan-06 


14 


1 


14 


1 


14 


1 


14 


1 




3C98 


0.030 


M 


Scp-06 


14 


1 


14 


1 


14 


1 


14 


1 




DA240 


0.036 
























3C192 


0.060 


S 


Nov-05 


60 


2 


60 


2 


120 


2 


120 


1 




4C73.08 


0.058 
























3C236 


0.101 


M 


Dcc-05 


14 


1 


14 


1 


14 


1 


14 


1 




3C277.3 


0.085 




Jan-06 


14 


1 


14 


1 


14 


1 


14 


1 




3C285 


0.079 


M 


Jan-06 


14 


1 


14 


1 


14 


1 


14 


1 




3C293 


0.045 


S 


Jan-06 


240 


2 


240 


3 


120 


3 


120 


2 


1.2 


3C305 


0.042 


M 


Apr-06 


14 


1 


14 


1 


14 


1 


14 


1 




3C321 


0.096 


S 
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14 
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14 


2 


30 
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30 


1 




3C326 


0.090 


S 


Mar-05 


60 


2 


60 


2 
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2 
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1 




3C382 


0.058 


S 


Aug-05 


60 


2 


60 


2 


120 


2 


120 


1 


1.2 


3C388 


0.092 


S 


Jul-05 


60 


2 


60 


2 


120 


2 


120 


1 


1.5 


3C390.3 


0.056 


S 


Aug-04 


240 


1 


240 


1 


120 


4 


120 


4 


1.0 


3C452 


0.081 


S 


Dcc-04 


60 


2 


60 


2 


120 


2 


120 


1 





Note. — Table presenting the basic parameters for the 3CRR sample as well as details of the Spitzer/IRS observations. 
Note that 3C403 and 3C445 overlap between the two samples, see Tabled for observation details of these two objects. 
Definitions are the same as Table [TI 



be varied, however, these three objects were all 
part of the same program (program 20719; P.I. 
Baum) and the observational patterns were iden- 
tical. All pointings had an integration time of 14 
sees, with 15 and 5 pointings for the SL and LL 
modules respectively. The middle number (i.e. 
8 in SL, and 3 in LL) pointing was centered on 
the target and, since these objects are relatively 
bright at MFIR wavelengths, this single central 
observation was used to extract the spectrum of 
each object. Because of the relatively short expo- 
sure times of the single exposures centered on the 
targets, these mapping mode observations tend 
to have lower S/N than the staring mode obser- 
vations used for the majority of the sample. In 
order to subtract the background, a median im- 
age of the other pointed observations was made, 
and we took this median image to represent the 
sky spectrum, subtracting this from the central 
pointing observation. To reduce the possibility of 
removing flux from the source in this process, we 
avoided using those pointed observations directly 
either side of the central pointing in the sky im- 
age (i.e. spectra 7 and 9 for SL, 2 and 4 for LL). 
Overall, the continuum level of the Mapping mode 
data matched well with the continuum measure- 
ment from Spitzer MIPS data measured over the 
24/im spectral filter response range (see below for 



further discussion). 

All the data were downloaded in the ba- 
sic calibrated data (BCD) format processed 
with the S18.7.0 data pipeline. These data 
were combined and subtracted using our own 
IDL code and cleaned using the SSC software 
IRSCLEAN MASK. The data were extracted us- 
ing SMART (v. 8. 1.2.) program dev eloped by the 



2004: Lebouteiller et al 



IRS Team at Cornell Univer sity (jHigdon et al 

2010[) . To extract the 



fluxes we employed the optimal extraction func- 
tion which uses a super sampled PSF and weights 
th e extracted spectra by the S/N of each pixel 



(|Lebouteiller et al.l [20100, assuming the objects 



are unresolved point sources at these wavelengths. 
The extracted spectra from both nod positions 
were compared by eye, and further cleaning by 
hand was performed where necessary i.e. by aver- 
aging over obvious strong sharp features that only 
appeared in one nod position. These spectra were 
then averaged together to produce the final spec- 
trum. Note that the IRS pipeline automatically 
accounts for variable slit losses with wavelength. 

17 of the 19 objects in the 3CRR sample have 
been observed with Spitzer/IRS. Two objects were 
not observed (4C73.08 and DA240) and the S/N 
of the data for 3C35 and 3C277.3 was too low to 
extract science quality data. Therefore 15/19 ob- 
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Fig. 1. — Spitzer/IRS spectra for the 2Jy sample. Common fine structure and H2 emission lines are indicated 
by vertical grey dashed lines. In objects with detected PAH features, the most prominent PAH features are 
indicated by blue line segments. Note the data for PKS0034-01 and PKS0035-02 were obtained from the 
Spitzer archive. Both potentially suffer from enhanced flux calibration uncertainties at short wavelengths, 
due to saturation of the peak- up detector (see Section [3] for details). 



jects (79%) have good IRS data. The data were re- 
duced in an identical manner to that described for 
the 2Jy sample; 5 objects were observed in Map- 
ping Mode and the rest in Staring Mode. Details 
of the observations can be found in Table [H 

For the majority of the resulting spectra (70%) 
the continuum level of the SL spectrum matched 
well with that of the LL spectrum. This demon- 
strates that the extraction techniques and back- 
ground subtraction method are robust for most 
objects, and that the sources are not spatially ex- 
tended. However, for 15 objects in the 2Jy sample 
and 2 objects in the 3CRR sample, the SL spec- 
trum and LL spectrum did not precisely match 



when we attempted to combine them together. 
The likely causes of this flux difference include the 
objects being partially extended, pointing errors, 
or poor calibration due to peak up saturation of 
the detector (discussed above). For most of these 
objects we applied a scaling factor to the SL spec- 
trum of between 0.8 and where the median 
scaling factor was 1.1 (see Tables [1] & [2]) . We kept 



■*For PKS0915-11 (Hydra A), PKS1938-15, 3C305, 3C382 
and 3C388 the scaling factor was higher, with values in the 
range 1.2 — 3.0. In these cases the objects may not have 
been exactly centered in the 3.6" SL slit because of pointing 
errors, or the mid-IR flux is significantly extended outside 
the SL slit. 
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Fig. 2. — Spitzer/IRS spectra for the 2Jy sample continued. Note the data for PKS0404+03 were taken in 
Mapping Mode and were obtained from the Spitzer archive. 



the flux level of the LL spectrum constant because 
the LL slit width is larger and therefore more likely 
to include the entire flux from any sources that are 
partially extended, or not perfectly centered. 

The spatially integrated 24/im flux values ob- 
tained from our Spitzer/MIPS photometric cam- 
paign can be used to check the flux calibration of 
the Spitzer/IRS data. The mean flux from the IRS 
spectra was calculated over the response range of 
the 24/xm filter for MIPS in the observed frame. 
We found that the estimated IRS fluxes agree to 
within 15% of the measured MIPS 24/im photo- 
metric values for all the sources, with a median of 
the ratio of to /i?S'24pm of 0.94. This 

provides evidence that our spectra capture the ma- 
jority of the mid-IR emission associated with any 
star formation in the host galaxies of the radio 



sources. 

Overall, we have useable IRS spectra for 56 
(89%) of the 63 objects in the combined 2Jy and 
3CRR sample. 

4. Data Analysis 

The fully reduced rest frame spectra for the 
43 observed/detected objects in the 2Jy sample 
and the 13 objects from the 3CRR sample (not 
in common with the 2Jy sample) are presented in 
Figures [1] through [HI The IRS spectral range is 5- 
38/im, however, at longer wavelengths the S/N is 
low. Therefore, we present the spectra only up to 
30/im in the observed frame and/or 27//m in the 
rest frame. 

The diversity in the spectra of the radio 
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Fig. 3.— Spitzer/IRS spectra for the 2Jy sample continued. Note the data for PKS0915-11, PKS0945+07 and 
PKS1602+01 were obtained from the Spitzer archive, and the data for PKS0947+07 were taken in Mapping 
Mode. PotentiaUy, the spectrum of PKS1602+01 suffers from enhanced flux calibration uncertainties at 
short wavelengths, due to saturation of the peak-up detector. 



sources shown in Figures [T] through |5] is strik- 
ing. Many of the objects show the characteristic 
features of AGN at mid-IR wavelengths, includ- 
ing strong fine structure emission lines, as well 
as weak si l icate emission or absorption features 
(|Hao et al.l 120071) . We have identified some of 



the most prominent fine-structure lines on the 
spectral plots: [Ne VI]A 7.65/im, [ArIII]A8.99/^m, 
[SIV]A10.5lAim, [NeII]A12.8lAim, [Ne V]A14.32^m, 
[NeIII]A15.56/im, [SIII]A18.7lAfm, [Ne V] A24.31//m, 
[OIV] A25.89/ini. We also indicate the location of 
H2 lines (S(l) through S(6)jl. It is noteworthy 



^Note [Aril] A6.99/im is not indicated in the plots and may 
be blended with H2 S(5)A6.91/^m. 



that nearly all the spectra show the high ionization 
[OIV]A25.89/im (i?io„ =54. 9eV) emission, when 
not redshifted outside the observable wavelength 
range. This indicates a high ionization state, as 
expected given the presence of powerful AGN in 
many of the sources. Low ionization potential 
lines such as [Nell] A 12.81 //m (£'io„=21.6eV) are 
also detected in many spectra. The contribution 
of starlight from the host galaxies is significant at 
short wavelengths (A < 8/im) in some low red- 
shift objects, seen as a sharp upturn in fiux at 
the blue end of the spectra. This is particularly 
apparent for objects with low power AGN, for ex- 
ample, the WLRG PKS1839-48 and PKS1954-55 
(see Figure SI). The fine structure, H2, and silicate 
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Fig. 4. — Spitzer/IRS spectra for the 2Jy sample continued. Note that the data for PKS1949+02 (Mapping 
Mode) and PKS2221-02 were obtained from the Spitzer archive. 



emission/absorption features wiU be discussed in 
depth in Paper II. 

A number of PAH dust emission bands exist 
within the spectral range of Spitzer/IRS («25), 
but at low spectral resolution they can be difficult 
to detect, for example the 12.7/im feature which 
blends with the strong [Nell] line at 12.81//m. 
Therefore, in this paper we focus on the three 
strongest PAH bands in the usable ZRS" wavelength 
range: the 6.2, 7.7, 11.3/Ltm features. The latter 
two PAH band are in fact blends of PAH emis- 
sion (7.60 & 7.85/xm, 11.23 & 11.33/im). Other 
bands that can make a significant contribution are 
the 8.6/im band and a blend of PAH emission at 
17/J,m (17.38 & 17.87^m), which can be strongly 
contaminated by H2 emission. 

To fit the data and measure the PAH fluxes 



we used PAHFIT vl.^, which is an IDL pro- 
gram developed by J.D.T Smith and B.T. Draine 
for studying PAH features in the mid-IR spectra 
of t he inner reg i ons o f local star-forming galax- 



ies ( Smith et all 120071 ) . The PAHFIT model is 



made up of 5 components: (1) starlight contin- 
uum, represented as a blackbody with temper- 
ature T=5,000K, (2) a featureless thermal dust 
continuum, represented by an array of 8 different 
black body continuum components with temper- 
atures from 35 to 300K, (3) pure rotational lines 
of H2, (4) fine structure forbidden emission lines, 
and (5) PAH emission features. 

We experimented with, and adapted, the PAH- 



S PAHFIT is made available under the terms of the GNU 
General Public License. 
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Fig. 5. — Spitzer/IRS spectra for the 2Jy and 3CRR sample objects. All the 3CRR IRS data were obtained 
from the Spitzer archive. The data for 3C98, 3C236, 3C277.3, 3C285 and 3C305 were taken in Mapping 
Mode and were obtained from the Spitzer archive. 



FIT code to check that it functioned appropri- 
ately for the samples of powerful radio-loud AGN 
presented in this paper. For example, we exper- 
imented with adding dust components into the 
model with black body temperatures 400K, 600K, 
lOOOK and 2000K, in order to account for the po- 
tential hotter dust continuum features in the ra- 
dio galaxy samples spectra. However, PAHFIT 
fits the data just as well without these extra hot 
components, which were therefore excluded in the 
final fit. In addition, when fitting the lO/im and 
18/im silicate absorption features, we found that 
the default PAHFIT extinction model did not fit 
the features well. After some experimentation 
it was found that the depth of the silicate ab- 
sorption feature was better fitted in models as- 



suming high ratios of 10 to 18/xm opacity. This 
implies that the 18/im feature is weak for our 
samples of powerful radio-loud AGN. Finally, the 
original PAHFIT model does not include silicate 
emission because galaxies with such emission were 
not included in th e original PAHFIT sample of 
Smith et al.l (j2007l) . However, for the 2Jy sam- 
ple, lO^m silicate emission is clearly detected at 
varying strengths in 19% of the observed/detected 
objects in the 2Jy and 3CRR samples. Therefore 
we made a further adaptation of the PAHFIT pro- 
gram to fit silicate emission as we ll as absorption 



follow ing the method employed in lGallimore et al 
()2010l) . 
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5. Results 



In this section we consider the rate of detection 
of PAH features in radio galaxies, and compare 
the various star formation diagnostics. 




Wavelength ^m) 



Fig. 6. — Spitzer/IRS spectra of the 3CRR sample 
continued. 



5.1. PAH detection 

It is important to distinguish PAH emission 
features from fine-structure emission lines, in or- 
der to make secure identifications of the PAH 
bands. Fortunately, the PAH emission bands ap- 
pear broader than the fine-structure lines, even 
in low resolution IRS spectra. The 11.3/im PAH 
feature falls in a region of the spectrum without 
strong fine-structure emission lines, making it an 
ideal choice for identifying the presence of PAH 
emission in the radio galaxy samples. In addi - 
tion, the study o f iDiamond-Stanic fc Riek3 (|2010[ ) 
found that the 11.3/zm feature is a more robust in- 
dicator of RSFA than the shorter wavelength PAH 
features in AGN host galaxies. Although it is more 
susceptible to the effects of silicate absorption and 
emission than shorter wavelength PAH features 
such as the 6.2/im feature, PAHFIT accounts for 
the effect of silicates when fitting the 11.3/im PAH 
emission. 

For the objects with good IRS spectra in the 
combined 2Jy and 3CRR sample, 9 (16%) ob- 
jects have PAH features that are clearly de- 
tected at high equivalent width (PKS0023-26, 
PKS0347-I-05, PKS0915-11, PKS1733-56, PKS2135- 
20, PKS2314 +03, 3C285, 3C293, 3C305). To 
show the PAH emission more clearly in these 9 
objects we present enlargements of their spectra 
between 5 and 12/im in Figures [7] and [8] for the 
2Jy and 3CRR samples respectively. 

Using the 11.3/im feature as a clean identifier, 
we find 8 more objects with low equivalent width 
(EW) or lower S/N PAH detections (PKS0806- 
10, PKS1151-34, PKS1814-63, PKS1934-63, 3C33, 
3C236, 3C321, 3C326). The detection strategy 
used to identify PAH emission in these objects 
was based first on checking that the candidate fea- 
ture was detected in both nod positions of the IRS 
spectra (where both were available). Second, we 
modeled each spectrum with PAHFIT, including 
or excluding the PAH emission components. Af- 
ter subtracting the PAHFIT model that did not 
include the PAH components from the spectrum, 
we were able to analyze the resulting residual spec- 
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Fig. 8. — : Magnified IRS spectra continued. Pre- 
senting the 3 objects in the 3CRR sample with vi- 
suahy identifiable PAH emission are shown. Note 
that 3C285 and 3C305 were observed in Mapping 
Mode and therefore suffer from lower S/N. 
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Fig. 7. — : Magnified IRS spectra focussing on the 
5-12 fj,Ta band of IRS to emphasize the dominant 
PAH bands at 6.2, 7.7 and 11.3/im. The 6 objects 
in the 2Jy sample with visually identifiable PAH 
emission are shown. 
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trum for evidence of PAH emission. 

To emphasize these lower level detections we 
present enlargements of the spectra in Figure [9] 
Below each spectrum in Figure IH] we plot the resid- 
uals from the PAHFIT model, where the solid 
line represents the residuals from subtracting the 
spectrum from the PAHFIT model that did not 
include PAH in the fit, and the dotted line rep- 
resents residuals of the PAHFIT model that in- 
cluded the PAH features. Four objects (PKS1814- 
63, PKS1934-61, 3C33, 3C326) have a prominent 
n.3^j.m PAH feature fitted weU by the PAHFIT 
model. In addition, low EW PAH features at 
7.7/im appear to be detected. Therefore we ar- 
gue that there is good evidence in these 4 objects 
for PAH emission. 

A further four objects (PKS0806-10, PKS1151- 
34, 3C236, 3C321) have lower EW or irregular 
PAH detections at 11.3//m. Because the 11.3/im 
feature also lies at the red end of the 10/im silicate 
absorption feature, potentially the wing of this fea- 
ture could be mistaken as a low EW signature of 
PAH. However, for PKS0806-10 and 3C321 the in- 
flection at 11.3/im is much sharper than expected 
for the wing of the 10/im silicate feature; in both 
objects the feature remains following subtraction 
of the PAHFIT model that includes silicate ab- 
sorption but excludes PAH features. Therefore we 
regard the detection of PAH features in these two 
objects as secure. 

PKS1151-34 also appears to have a prominent 
11.3/Ltm feature, albeit with an irregular shape. 
The irregular shape is due to the fact that the fea- 
ture coincides with the join between the SL and 
LL parts of the spectrum which hinders a more 
definitive detection. 

Finally we note that 3C236 also shows a pos- 
sible 11.3/im feature detected at low S/N in its 
low resolution IRS spectrum (see Figure 9). Al- 
though this detection in not secure based on the 
low resolution IRS spectrum alone (taken in map- 
ping mode), it has been confirmed with a much 
higher S/N spectra taken using the high resolution 
mode of Spitzer/IRS (P. Guillard, private commu- 
nication) . 

5.2. PAH emission strength 

PAHFIT uses a Drude profile method to recover 
the full strength of the PAH dust emission features 



and blends. Drude profiles are known to recover 
the flux in PAH features more accurately than 
methods which estimate the underlying continuum 
using line segments or spline curve fits, because 
they better account for the power in the wings 
of the broad emission profiles. PAHFIT also cor- 
rects the PAH fiuxes for the corresponding effects 
of the silicate absorption model fitted to the data. 
Converting these fluxes into luminosities gives an 
opportunity to compare the relative strengths of 
the emission features for all the objects with IRS 
spectra. 

Our PAHFIT model fits the detected 11.3//m 
PAH emission features and returns flux values for 
all the objects in the 2Jy and 3CRR samples with 
PAH detections. However, in the case of those 
objects in which the PAH features are not de- 
tected, it is necessary to derive upper limits. To 
do this in a robust fashion, a scaled IRS spec- 
trum of PKS0915-11 (Hydra A), which is domi- 
nated by prominent PAH features, was added to 
the spectrum of each object without evidence for 
PAH emission. The scaling factor was varied until 
a PAH feature at 11.3/im was just detected in a vi- 
sual inspection of the combined spectrum. Again 
the 11.3/im feature was used, as this is not con- 
taminated by flne structure lines. Multiplying the 
measured PAH flux for PKS0915-11 by the scaling 
factor then gives a robust upper limit on the flux 
of the PAH. The PAH fiuxes and upper limits are 
presented in Tables [3] and |31 

Using these upper limits we can test whether 
the objects with detected PAH emission have a 
unique signature for starbursts, in that the lumi- 
nosity of the detections is greater than the upper 
limits, or whether significant PAH emission could 
exist in all the spectra but may not be detected 
due to low S/N data. In Figure ITOl we plot the 
11.3/im luminosities against the AGN power indi- 
cator i[oiii]- From this figure it is clear that, on 
average, the PAH luminosities of objects with high 
EW PAH detections are higher than the majority 
of PAH upper limits for a comparable i[o///] 
AGN power. 

The difference in Figure [10] between the PAH 
detections and the upper limits is clearest for low 
power AGN, but becomes less clear for higher 
power AGN. Two factors may contribute to this 
trend. First, because of limits on the exposure 
times, some of the more powerful, distant ob- 
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Fig. 9. — : Magnified IRS spectra focussing on the 5-12/ini band presenting the objects deemed to have low 
EW PAH. Below each spectrum are plotted the residuals from the PAHFIT fit , the solid line representing 
the PAHFIT fit without PAH and the dotted line representing a fit with PAH component. See Section [5TT] 
for more details. 
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jects have lower S/N spectra than their lower 
power counterparts at low redshifts, leading to 
higher upper limit flux values. In this context it 
is notable that the group of upper limits around 
LpAH ~ 10^^ W, L[oiii] 10^^ W — which in- 
cludes PKS0409-75, PKS1306-09, and PKS0252- 
72 — includes some of the highest redshift objects 
in the 2Jy sample (z > 0.5); these objects have 
relatively low S/N IRS spectra. Second, the rela- 
tively stronger AGN dust continuum emission in 
the more powerful objects could also lead to higher 
upper limits for the PAH emission. 

The objects with low EW PAH detections are 
also marked in Figure 1101 In general, they occupy 
the same region of the plot as the upper limits, 
but tend to lower PAH luminosities compared to 
the higher EW PAH detections. 

From these data alone it is impossible to en- 
tirely rule out RSFA activity, at some level, in the 
objects with PAH upper limits. However, the posi- 
tion of the upper limits for lower luminosity AGN 
argues against powerful starburst activity in these 
objects. 
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Fig. 10.— : Plot of LpAHii.3 versus i[oiii] 

By combining the spectra that have no evidence 
for PAH features, it is possible to create a spec- 
trum with high signal-to-noise to investigate the 
possibility of low level PAH emission in the ob- 
jects without individual detections. We present 
median spectra for the objects without PAH de- 
tections in the combined 2Jy and 3CRR sample 
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Fig. 11. — : Median spectra for the combined 2Jy 
and 3CRR samples; objects with PAH detections 
(blue), objects without PAH detections (black), 
BLRG without PAH detections (red), NLRG with- 
out PAH detections (green). The median spectra 
were normalized to IJy at 20/im and offset in the 
plot in order to compare them. 

in Figure [TT] These median spectra were created 
by re-sampling the spectral data in wavelength 
bins and normalizing all the spectra to have the 
same flux at 20/im. While the median spectrum 
of the PAH identified objects shows the character- 
istic PAH emission dominating the spectrum, the 
median spectrum for objects in which we have not 
individually identified PAH shows no evidence for 
PAH features. 

Also in Figure [TT] we present median spectra 
for objects without individually identified PAH 
features divided into BLRG and NLRG. We can 
use these median spectra to test the idea that a 
stronger AGN continuum in BLRG may mask low 
level PAH emission, which in turn may be more 
easily detected in NLRG. The test is again neg- 
ative, revealing no evidence for PAH emission in 
the NLRG. 
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Table 3 

The 2Jy Sample - PAH Flux Data 



PKS Name 


other 




z 


Opt. 


Rad. 


Sii.3(W/m^ 


) 


70/24 


UUZo — 




0. 


322 






l^o.4 It U. X 


10 


■17 


12.2 i 2.1 


UU04 — Ul 




n 
u. 


U ( o 


VV ijxt,^ 


FRII 


< 6.Dr!j-lo 








uuoo — 




n 
u. 


220 


BLR.G 


ri.ii ) 


< 1.3E-17 






1 Q -t- n /I 
i.y m U.4 


uuoo-i-uy 




n 
u 


.loo 


BLRG 


FRII 


< O.Orj-lo 






1 2 ih 2 


uuoy — 44 




n 
u. 


.04D 


NLRG 


FRII 


^ o nT? 1 7 






9 1 -i- n 1 
z. 1 m u. 1 


UU40 — 4^ 







116 


WT Rr* 


FRII 


^ f< fiT? 1 S 
<^ D.DIi;-10 






9 + 03 


UlUO — ID 




n 

u. 


A nn 

.4UU 


NLR,G 


FRII 


< 2.0E-17 






^ 1 


0117—15 


3C38 


0. 


.565 


NLRG 


FRII 


< 6.6E-18 






3.3 lb 0.3 


0213 — 13 


3C62 





.147 


NLRG 


FRII 


< 6.6E-18 






0.9 ± 0.1 


0235-19 


OD-159 


0. 


.620 


BLRG 


FRII 


< 1.3E-17 






1.3 ± 0.2 


0252-71 




0, 


.566 


NLRG 


CSS 


< 6.6E-18 






< 3 


0347+05 




0. 


.339 


WLRG 


FRII 


(8.9 ± 0.1) X 


10" 


17 


8.8 ± 1.3 


0349—27 







.066 


NLRG 


FRII 


< 6.6E-18 






4.8 ± 0.2 


0404+03 


3C105 


0. 


.089 


NLRG 


FRII 


< 2.0E-17 






2.3 ± 0.1 


0409-75 




0, 


.693 


NLRG 


FRII 


< 9.9E-18 






7.3 ± 1.9 


0442 — 28 




0. 


.147 


NLRG 


FRII 


< 6.6E-18 






1.4 ± 0.2 


0620-52 




0, 


.051 


WLRG 


FRI 


< 6.6E-18 






10.5 ± 0.4 


0625-35 


OH-342 


0, 


.055 


WLRG/BLLac 


FRI 


< 6.6E-18 






1.8 ± 0.1 


0625-53 




0, 


.054 


WLRG 


FRII 


< 6.6E-18 






< 6 


0806-10 


3C195 


0. 


.110 


NLRG 


FRII 


(1.8 ± 1.5) X 


10" 


17 


1.9 ± 0.01 


0859-25 




0, 


.305 


NLRG 


FRII 


< 6.6E-18 






0.9 ± 0.3 


0915-11 


Hydra A 





.054 


WLRG 


FRI 


(1.1 ± 0.1) X 


10" 


16 


12.9 ± 0.6 


0945+07 


3C227 


0, 


.086 


BLRG 


FRII 


< 4.0E-17 






0.4 ± 0.1 


1136-13 




0, 


.554 


Q 


FRII 


< 2.0E-17 






1.7 ± 0.2 


1151-34 




0. 


.258 


Q 


CSS 


(3.0 ± 0.1) X 


10" 


17 


3.2 ± 0.2 


1306-09 




0, 


.464 


NLRG 


CSS 


< 6.6E-18 






4.7 ± 0.5 


1355-41 




0. 


.313 


Q 


FRII 


< 1.3E-17 






1.2 ± 0.04 


1547-79 




0, 


.483 


BLRG 


FRII 


< 1.3E-17 






2.4 ± 0.2 


1559+02 


3C327 


0, 


.104 


NLRG 


FRII 


< 2.6E-17 






1.9 ± 0.02 


1602+01 


3C327.1 


0, 


.462 


BLRG 


FRII 


< 1.3E-17 






1.6 ± 0.3 


1648+05 


Here A 


0, 


.154 


WLRG 


FRI 








< 9 


1733-56 




0, 


.098 


BLRG 


FRII 


(1.3 ± 0.1) X 


10" 


16 


5.2 ± 0.1 


1814-63 




0. 


.063 


NLRG 


CSS 


(8.2 ± 0.5) X 


10" 


17 


2.3 ± 0.0 


1839-48 




0, 


.112 


WLRG 


FRI 


< 6.6E-18 






3.5 ± 0.9 


1932-46 




0. 


.231 


BLRG 


FRII 


< 6.6E-18 






7.1 ± 0.7 


1934-63 




0, 


.183 


NLRG 


GPS 


(1.3 ± 0.1) X 


10" 


17 


1.1 ± 0.1 


1938-15 




0, 


.452 


BLRG 


FRII 


< 1.3E-17 






2.9 ± 0.6 


1949+02 


3C403 


0, 


.059 


NLRG 


FRII 


< 4.0E-17 






1.8 ± 0.02 


1954-55 




0, 


.060 


WLRG 


FRI 


< 6.6E-18 






3.3 ± 1.1 


2135-14 




0, 


.200 


Q 


FRII 








1.1 ± 0.05 


2135-20 


OX-258 


0, 


.635 


BLRG 


CSS 


(1.1 ± 0.1) X 


10" 


16 


8.7 ± 1.0 


2211-17 


3C444 


0, 


.153 


WLRG 


FRII 








< 18 


2221-02 


3C445 


0, 


.057 


BLRG 


FRII 


< 6.6E-17 






0.8 ± 0.02 


2250-41 




0, 


.310 


NLRG 


FRII 


< 2.6E-18 






1.9 ± 0.2 


2314+03 


3C459 


0, 


.220 


NLRG 


FRII 


(1.1 ± 0.1) X 


10" 


16 


10.3 ± 0.1 


2356-61 




0, 


.096 


NLRG 


FRII 


< 1.3E-17 






1.8 ± 0.1 



Note. — Tabic presenting luminosities derived from the fluxes extracted using the PAHFIT 
program for all 2Jy objects with useful IRS data. Uncertainties are derived from the PAHFIT 
model fit. See Tabl({6lfor ratios of 11.3/6.2^m and 11.3/7.7Mm fluxes. PKS0347+05 had been 
classified in our previous work as a BLRG (TOT, DOS, D09), however, we have since deter- 
mined that the broad-line component comes from a companion galaxy and not the radio galaxy, 
therefore it is now classified as a WLRG. 
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Table 4 

The 3CRR Sample - PAH Flux Data 



Name 


z 


Opt. 


Rad. 


S'll.3(W/ 




70/24 




U.UDU 


NLRG 


FRII 


l^o. y i u.Dj X 




1 -1- n m 
1.0 Hz u.uo 


3C35 


0.067 


WLRG 


FRII 






21.0 ± 8.6 


3C98 


0.030 


NLRG 


FRII 


< 1.3E-17 




0.8 ± 0.1 


DA240 


0.036 


WLRG 


FRII 






8.3 ± 1.5 


3C192 


0.060 


NLRG 


FRII 


< 2.6E-17 




2.4 ± 1.1 


4C73.08 


0.058 


NLRG 


FRII 






0.5 ± 0.1 


3C236 


0.101 


WLRG 


FRII 


(5.4 ± 0.1) X 


10-1^ 


3.7 ± 0.3 


3C277.3 


0.085 


WLRG FRI/FRII 






2.1 ± 0.4 


3C285 


0.079 


NLRG 


FRII 


(1.3 ± 0.1) X 




4.3 ± 0.1 


3C293 


0.045 


WLRG FRI/FRII 


(2.3 ± 0.1) X 




9.7 ± 0.2 


3C305 


0.042 


NLRG FRII/CSS 


(1.4 ±0.3) X 


10-1'^ 


7.1 ± 0.1 


3C321 


0.096 


NLRG 


FRII 


(1.5 ± 0.1) X 


10-1^ 


3.4 ± 0.02 


3C326 


0.090 


NLRG 


FRII 


(1.1 ± 0.1) X 


10-1^ 


< 13 


3C382 


0.058 


BLRG 


FRII 


< 1.3E-17 




0.6 ± 0.04 


3C388 


0.092 


WLRG 


FRII 


< 6.6E-18 




< 4 


3C390.3 


0.056 


BLRG 


FRII 


< 2.6E-17 




0.8 ± 0.01 


3C452 


0.081 


NLRG 


FRII 


< 6.6E-18 




1.0 ± 0.1 



Note. — Tabic presenting luminosities derived from the fluxes extracted 
using the PAHFIT program for all 3CRR the objects with useful /i?S' data. 
Note that 3C403 and 3C445 overlap between the two samples; sec Tabic m 
for PAH data for these latter objects. 



5.3. Comparison with optical star forma- 
tion indicators 

Careful spectral synthes is modeling of deep op- 



tical spectr a for the 2Jy ('Ta dhunter et al 



3CRR (see Sicken et al J (|20inr and references 
therein) samples has allowed identification of ob- 
jects with evidence for RSFA at optical wave- 
lengths. 

Considering the objects with high EW PAH 
detections that also have deep optical spectra 
available, it is notable that all 7 of these ob- 
jects show evidence of young stellar populations in 
their optical spectra. In addition, we have iden- 
tified RSFA through PAH detection in 2 objects 
(PKS0347-H05, PKS1733-56) for which the con- 
tamination by direct AGN emission components 
precluded the identification of any young stellar 
populations at optical wavelengths. 

Out of the 12 objects with good optical evi- 
dence for RSFA and IRS spectra, only 3 objects 
show no evidence for PAH features (PKS0409-75, 
PKS0620-52, PKS1932-46). The relatively low 
S/N of the spectra in the region of the 11.3/im 
PAH feature for these 3 objects (S/N = 2, 13 and 
4 respectively, compared with the sample average 



Wills et all ^04, 2008; 



2002 



20071 ) and 



of 24) could help to explain why no PAH emission 
is detected. 

Overall the optical and PAH evidence for RSFA 
correlates well. The correlation is unlikely to be 
perfect since the two star formation indicators do 
not necessarily sample star formation activity on 
the same spatial scales, or at the same evolution- 
ary stages of the starbursts. For example, in a 
significant subset of the radio galaxies with opti- 
cally identified young stellar populations we are 
detecting post-starburst stellar populations that 
are observed a significant period (~ 0.1 — 2 Gyr) 
after t he peak of merger induced sta r formation ac- 
tivity (|Tadhunter et al.ll2005l boill ). In addition, 
identification of RSFA at optical wavelengths is 
not possible in some objects due to a strong AGN 
continuum that masks these signatures. Equally, 
PAH identification may be hindered by low S/N 
spectra and/or calibration errors such as noted for 
3C15 (see Section [3]). 

Combining the objects that have high EW PAH 
detections (9) and the objects that have clear op- 
tically identified RSFA (12) makes a total of 14 
objects (22%) with strong evidence for RSFA in 
the combined 2Jy and 3CRR sample. 

Next we consider the 6 objects with low EW 
PAH detections that also have deep optical spectra 
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available. Only 2 of these objects - 3C236, 3C321 
- have evidence for RSFA at optical wavelengths. 
The low EW PAH detections in the remaining 
4 objects - PKS0806-10, PKS1934-63, 3C33 and 
3C326 - may represent a low level of star forma- 
tion that could not be detected at optical wave- 
lengths. However, it is also important to consider 
that the low EW PAH emission may originate from 
ISM diffuse PA H emission rather th an be associ- 
ated with RSFA (jKaneda et al.ll2008l ). particularly 
for objects with low AGN continuum strength such 
as 3C326. 

Finally, we return to the question posed in the 
introduction of whether there is a substantial pop- 
ulation of radio galaxies in which the RSFA is 
missed at optical wavelengths due to dust obscu- 
ration. We find that, of the 36 objects that have 
good optical spectroscopic information and show 
no sign of RSFA in their optical spectra, only the 
4 objects (11%) discussed above - PKS0806-10, 
PKS1934-63, 3C33 and 3C326 - have PAH fea- 
tures detected in their mid-IR spectra; and these 
are all low EW detections. Therefore, the inci- 
dence of objects in our samples that have obscured 
RSFA that is not already detected at optical wave- 
lengths is relatively minor. 

5.4. Mid- to far-IR color 

A common method for identifying star formation 
activity involves using the mid- to far-IR (MFIR) 
color. Both the 2Jy and 3CRR samples have com- 
plete Spitzer MIPS data which can be used to in- 
vestigate the F(70)/F(24) flux ratio (D09, DIO). 
Empirically, starburst galaxies are known to be as- 
sociated with cool MFIR colors, as expected in the 
case of illumination of extended dust structures by 
spatially distributed star forming regions. Unlike 
PAH emission, the MFIR color diagnostic is not a 
clean identifier of RSFA, because the far-IR emit- 
ting dust can be heated by AGN and/or young 
stellar populations. However, in our previous work 
(D09, DIO) we have shown that optically identified 
RSFA objects do tend to have cooler MFIR col- 
ors. With the addition of the more complete PAH 
identified RSFA objects from the 2Jy and 3CRR 
samples, we can thoroughly test the MFIR color 
star formation diagnostic for radio-loud AGN. 

Figure [12] presents a histogram of the MFIR 
colors for the 2Jy and 3CRR samples, with the 
RSFA objects indicated. It is immediately clear 



that all but one of the 11 objects with F(70) /F(24) 
>5 show evidence for RSFA on the basis of opti- 
cal and/or PAH data, whereas the majority of ob- 
jects lacking such evidence cluster around MFIR 
colors of F(70)/F(24) « 0.5-2.5. However, there is 
a mix of objects in the range 2.5<F(70)/F(24)<5, 
which may represent a MFIR color transition re- 
gion between AGN- and RSFA-dominated objects. 
Objects in this transition region with F(70)/F(24) 
^4.5 include PKS0349-27 and PKS1306-09. We 
note that these latter two objects do not have the 
high quality spectra required to detect the signa- 
tures of RSFA at optical wavelengths. It is also 
noteworthy that PKS1306-09 is a CSS object that 
shows a strong far-IR excess, in the sense that it 
falls well above the L7o^m vs i[oiii] correlation 



(see Section WM . 

All the objects with either high EW PAH 
and/or optical evidence for RSFA (highlighted in 
dark blue in Figure [T^ show relatively cool MFIR 
colors with 70/im/24/im > 3. Therefore this study 
finds that MFIR color is an excellent indicator 
of starburst activity in radio galaxies^. This is 



consis tent with the results found by ,Brandl et al 



(120061) for AGN u s ing th e F(15)/F(30) ratio, and 
bv IVeilleux et al.l (l20Q9l) using the F(60)/F(25) 
ratio for at large sample of ULIRGs with varying 
AGN contribution. 

Of the 8 objects with low EW PAH detections, 
5 have been identified in Figure I12[ while a fur- 
ther two objects (3C236 and 3C321) are already 
highlighted in this figure as optical RSFA objects 
(dark blue) with F(70)/F(24) = 3.7 and 3.4 respec- 
tively. The eighth low EW PAH object (3C326) is 
not included in this figure as it has an upper limit 
on F(70)/F(24). The figure shows that the MFIR 
colors of the low EW PAH objects are similar to 
those of objects lacking PAH detections. This re- 
sult is not surprising assuming that the low EW 
PAH detections represent a low level of star for- 
mation activity or indeed if the PAH has a diffuse 
ISM emission origin rather than a stellar origin 
as mentioned in above in Section [5.3l It is notable 
that the two low EW PAH objects with the coolest 
colors, 3C236 and 3C321, mentioned above, also 
have optical RSFA detections. 

It is also important to consider the possibility 
of RSFA in the 3 objects from the 2Jy sample and 
the 4 objects from the 3CRR sample that were 
not observed/detected with Spitzer/IRS. None 
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70|im/ 24|im 



] no-RSFA objects 
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] low Eqw PAH objects 



■ II II 



Fig. 12. — : Histogram of 70/xm/24/im color for the combined 2Jy and 3CRR sample. Note that the figure 
does not include the 5 objects in the 2Jy sample and the 2 objects in the 3CRR sample with upper limits 
on their 70/im flux. 



of these objects have optical evidence for RSFA. 
However, two of the closest objects (z«0.03) 3C35 
and DA240 have cool MFIR colors (=21.0±8.6 
and 8.3±1.5 respectively). We have excluded 
3C35 from this analysis because its measured 
F(70)/F(24) has a large uncertainty. DA240 is the 
only object with cool MFIR colors F(70)/F(24)>5 
without clear optical or PAH identification of 
RSFA. On the other hand, this object does show 
evidence for RSFA in the form of a far-IR excess 
(see Section [53] and Table [S]). 

5.5. Far-IR excess 

Our previous work has shown that the objects with 
optically identified RSFA have enhanced emission 
at far-IR wavelengths compared to the majority 
of the sample objects (T07, D09, DIO). Several 
objects in these previous studies also showed ap- 
parent enhancements in their far-IR emission, but 
optical identification of RSFA was not possible. 
We are now able to revisit this investigation with 
the benefit of the IRS data and PAH detections to 
see if our previous result holds with the new, more 
complete, RSFA identifications. 

In Figure [T2] we plot L24/jm and £70^™ against 
-^[Oiii] foi' 2Jy and 3CRR samples. As dis- 
cussed in Section [TJ L[oiii] is a good indicator 
of the intrinsic power of the AGN. Consequently, 
we interpret the similar slopes of the correlationf|j 
between L24^m, iro^m and L[oiu] seen in Figure 



[T51 as direct evidence for AGN illumination as the 
dominant heating mechanism of the thermal mid- 
to far-IR emitting dust (D09, DIO). Considering 
the spectral diversity seen in the IRS data pre- 
sented in this paper, the correlation between the 
mid-IR (L24pm) and L[oiii] is striking. 

Following the method of our previous investi- 
gations we have marked with blue stars in Fig- 
ure [13] those objects with high EW PAH detec- 
tions and/or optical evidence for RSFA in the 2Jy 
and 3CRR samples. Consistent with our previ- 
ous studies, there is little difference between the 
distributions of RSFA and non-RSFA objects in 
the L2i^m vs i[oiii] plot (Figure [13] (a)). In con- 
trast. Figure [TS] (b) clearly shows that all but 
one (PKS1932-46jl of the 14 objects with opti- 
cal and/or PAH evidence for RSFA in the 2Jy and 
3CRR samples lie above the L^Q^rn versus i[oiii] 
correlation defined by the regression line fitted to 
the no-RSFA objects. 

This boosting of the objects with evidence for 



RSFA above the correlation between L 



70fi.m 



and 



^[Oiii] in Figure [13] — by more than an order 
of magnitude in some cases — provides strong 



^Spearman partial rank correlation tests presented in DIO 
show the correlations are intrinsic and do not arise because 



^[OIII] ^^'^ Lmfir are independently correlated with red- 
shift. The null hypothesis that the variables are unrelated 
is still rejected at the >99.5% level of significance. 
^This object is discussed in detail in DIO. Amongst the 
objects with optically identified RSFA, this object has the 
weakest evidence for RSFA in its near-nuclear regions, with 
the evidence confined to the detection of HII region-like 
emission line ratios in th e extended halo of the galaxy 
l IVillar-Martm et al.|[2005f) . 
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Fig. 13. — : Plot showing i24pm and LiOfj^m versus i[oiii] where those objects identified to contain evidence 
for RSFA with optical and/or PAH diagnostics are indicated with blue stars. The lines shown are the 
bisectors for no- RSFA objects of linear least squares fits oi x on y and y on x. The fit also does not include 
the 7 objects with upper limits in either [OIII] or 70/im. The 24^m, 70/xm and [OIII] fiuxes an d luminosities 
for th e 2Jy sample, can be found in DOS and D09; fluxes for the 3CRR sample are presented in lPicken et al.l 

(|2nin[ ). 
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evidence for enhancement of the thermal far-IR 
emission by RSFA. Conducting a one dimensional 
Kolmogorov-Smirnov two sample test on the com- 
bined sample, comparing the vertical displace- 
ments from our fitted regression line in the Ljo^rn 
vs i[oiii] plot, we find that we can reject the 
null hypothesis that the objects with and without 
independent optical/PAH evidence for RSFA are 
drawn from the same parent population at a bet- 
ter than 0.01% level of significance. Note that, in 
conducting this analysis, we included upper limits 



on L 



70fJ.m 



and L 



[OIII] 



as actual measurements. 



To take into account the upper limits in a more 
rigorous manner, we also investigated the differ- 



ences in the distributions of L 



70iim 



and L 



[OIII] 



including the upper limits, usi ng survival statis- 



tics in the ASURV package (jlsobe et al.l Il986t 
iLavallev et al. 1992 ) as implemented in IRAF. 
Specifically, the "twosampt" task, which com- 
putes several nonparametric two-sample tests for 
comparing two or more censored data sets, was 
used. These tests revealed that the 70/im flux 
was significantly different in the starburst ver- 
sus the non-starburst populations. For instance, 
the Gehan's Generalized Wilcoxon test estimated 
that the probability (p) that the two sets of cen- 
sored data belonged to the same population was 
= 0.023 (test statistic = 2.268) while the Logrank 
test estimates p = 0.013 (test statistic — 2.483). 
On the other hand, the [OIII] flux was not sig- 
nificantly different between the two populations: 
the Gehan's Generalized Wilcoxon test gave a test 
statistic — 0.531 and p — 0.595, while the Logrank 
test gave a test statistic ~ 0.123 and p — 0.902. 

We have also plotted the objects with low EW 
PAH detections on Figure[T31 Many of the low EW 
PAH objects fall above the Ljo^m versus i[oiii] 
correlation, but the result is not as clear-cut as it 
is for the high EW PAH objects. It is noteworthy 
that the low EW PAH object furthest above the 
fitted regression line in the L-jt^^m vs -Zj[oiii] corre- 
lation - 3C321 - has the coolest MFIR colors out 
of all the low EW PAH objects (F(70) /F(24) = 
3.4) and also optical identification of RSFA. This 
object is indicated with a filled star in Figure fT3l 

In terms of whether the far-IR excess alone is 
a good indicator of RSFA, we find that all 13 of 
the objects that fall more than 0.5 dex (factor of 
3) above the Lro/jm and L[oiii] regression line also 
show evidence for significant RSFA based on at 



least two of the other indicators. We conclude that 
far-IR excess is as useful an indicator of RSFA in 
radio galaxies as MFIR color. 

6. Discussion 

6.1. Star formation activity in radio-loud 
AGN 

Our multi-wavelength dataset allows us to es- 
timate the maximum percentage of objects that 
show evidence for RSFA in the combined 2Jy and 
3CRR sample by considering all the main optical 
and infrared diagnostics of star formation activ- 
ity. In order to summarize the results we have 
collated the data for the various RSFA identifiers 
in Table [5l Several objects are not complete in 
terms of the available optical, PAH, far-IR excess 
and color tests of RSFA (indicated by a "U" in the 
Table). However, all the objects have at least two 
separate pieces of RSFA diagnostic information. 
We further note that, of the 17 objects with opti- 
cal, low or high EW PAH, color-based or infrared 
excess identification of RSFA, 7 objects are iden- 
tified as RSFA based on all four diagnostics, and 
a further 9 objects on the basis of two diagnostics. 

Considering now the full 2Jy-|-3CRR sample 
of 63 objects: 19% (12 objects) show unambigu- 
ous spectroscopic evidence for RSFA at optical 
wavelengths; 21% (13) have cool MFIR colors 
(F(70/F(24) > 5); 14% (9) have high EW PAH 
emission; 22% (14) lie more than 0.5 dex (a fac- 
tor of 3) above the regression line fitted to the 
L7o^m versus L[oiii] correlation. Overall we find 
that 27% (17) objects show at least one of these 
indicators. Including in the analysis the objects 
with low EW PAH detections (5 additional objects 
not identified using the other diagnostics) , we esti- 
mate that the maximum percentage of radio-loud 
AGN showing evidence for RSFA in the combined 
2Jy and 3CRR sample is 35%. 

This combined result from the various diagnos- 
tic methods provides conclusive evidence that a 
one to one connection between RSFA and AGN 
activity does not exist for the majority of interme- 
diate redshift radio-loud AGN. This conclusion is 
consistent with our previous investigations, as well 



(Tadhunter et al. 20071 D08: D09: 


Clearv et al. 


20071: IShi et al.l 2007: Fu & StocktonI 


20091) 



In a study of 33 3C quasars and radio galaxies 
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Table 5 

Recent star formation activity detections in the 2Jy and 3CRR samples 



Name 


other 


Optical 




L^oior 


r ar-lrv 


Name 


other 


Optical 


IT An. 


k^oior 


rar-irv 


0023-26 




✓ 


✓ 


✓ 


✓ 


1814-63 




U 


low 




✓ 


0034-01 


3C15 










1839-48 












0035-02 


3017 


- 


- 


- 


- 


1932-46 




✓ 


- 


✓ 


- 


0038+09 


3C18 










1934-63 






low 






0039-44 












1938-15 












0043-42 












1949+02 


3C403 










0105-16 


3032 






u 




1954-55 












0117-15 


3038 










2135-14 




U 


U 






0213-13 


3062 










2135-20 


OX-258 


✓ 


✓ 


✓ 


✓ 


0235-19 


OD-159 










2211-17 


3C444 




U 


U 




0252-71 








u 




2221-02 


3C445 










0347+05 




U 


✓ 


✓ 


✓ 


2250-41 












0349-27 




U 








2314+03 


3C459 


✓ 


✓ 


✓ 


✓ 


0404+03 


3O105 










2356-61 












0409-75 




✓ 




✓ 


✓ 


3C33 






low 






0442-28 












3C35 






U 


u 




0620-52 




✓ 




✓ 




3C98 












0625-35 


OH-342 










DA240 






U 


✓ 


✓ 


0625-53 








u 




3C192 












0806-10 


30195 




low 






4C73.08 






U 






0859-25 












3C236 




✓ 


low 




✓ 


0915-11 


Hydra A 


✓ 


✓ 


✓ 


✓ 


3C277.3 






U 






0945+07 


30227 


u 








3C285 




✓ 


✓ 


✓ 


✓ 


1136-13 




u 








3C293 




✓ 


✓ 


✓ 


✓ 


1151-34 




u 


low 






3C305 




✓ 


✓ 


✓ 


✓ 


1306-09 




u 






✓ 


3C321 




✓ 


low 




✓ 


1355-41 




u 








3C326 






low 


u 




1547-79 




u 








30382 




u 








1559+02 


30327 










30388 








u 




1602+01 


30327.1 










3O390.3 




u 








1648+05 


Here A 




U 


u 




3C452 












1733-56 




u 


✓ 


✓ 

















Note. — Tabic showing the RSFA detections for the 2Jy and 3CRR samples. A dash means no evidence for 
RSFA activity exists for that diagnostic method. Optical column: a tick indicates that the object have young 
stellar populations identified in careful modeling of optical spectra, U (uncertain) identifies objects where it was 
not possible to model the spectrum effectively. PAH column: a tick indicates positive strong PAH detection, low 
indicates a low EW or low S/N PAH detection, while U indicates that good IRS spectra are not available. Color 
column: a tick indicates that an object has 70yiim/24/j,m>5, U indicates that the object has an upper limit in 70jum 
flux. Far-IR column: a tick indicates that the object has a far-IR excess in the sense that it falls 0.5 dex (a factor 
of 3) above the regression line fitted to LjQ^nx versus L^qh^ correlation, not including object with upper limits. 
Note that two of the objects in the 3CRR sample are in common with the 2Jy sample: 3C403 (PKS1949+02), 
3C445 (PKS2221-02). 
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at int ermediate redshifts (0.5<z<l.l) IClearv et al 
20071 concluded that star formation does not con- 
tribute significantly to the MFIR emission. This 
conclusion was based on the 15/im/3G/im mid-IR 
flux ratio also use d in other investigations (e.g. 
Brand! et ablbOOd ). 



The study of lShi et all (|20071 ) used Spitzer IRS 
data to investigate PAH emission in SCR radio 
galaxies, as well as other samples of AGN. Eight 
objects overlap with the 2Jy an d 3CRR samples 



presented here. The results of IShi et al.l (|2007t ) 



agree with our study, detecting PAH features in 
only 2 of the overlapping objects (3C293, 3C321). 
The exception is 3C33, where we have identified 
low EW P AH e mission that was not identified in 
(|2007l) . Moreover, only one o bject out 



Shi et al 



of the 94 3CR objects included in the IShi et al 



(|2007f) study has both 7.7 and 11.3/im PAH fea- 
tures detected. Although their sample covers a 
wider redshift range (z<1.5), and their PAH detec- 
tion method is likely less sensitive than our study, 
the results shows little evidence for powerful dom- 
inant RSFA components in 3CR objects. 

Finally, IFu fc StocktonI (|2009h studied a sample 
of 12 FRII radio-loud quasars with z«G.3, a rela- 
tively similar sample of objects in terms of redshift 
to the 2Jy sample presented in this paper. They 
found no evidence for PAH emission in any indi- 
vidual object, or from combined spectra. 

Clearly, based on several studies of the PAH 
features, as well as the other diagnostics discussed 
in this paper, there is little evidence for powerful 
RSFA activity in the majority of radio-loud AGN. 

We emphasize that, due to the radio flux lim- 
ited selection method for the two samples, the 
higher redshift objects are likely to be more pow- 
erful radio-loud AGN. It is conceivable that this 
bias could affect the detection rate of PAH or other 
RSFA diagnostics across the redshift range of the 
sample. However, we find no obvious trend with 
redshift in terms of the rate of detection of RSFA. 
For example, the rate of detection of RSFA ac- 
tivity in the low redshift 3CRR sample (z < 0.1: 
36% have optical and/or PAH detection evidence, 
including low EW PAH detections) is not signifi- 
cantly different from that in the intermediate red- 
shift 2Jy sample (0.05 < z < 0.7: 28%). This 
confirms that the diagnostic methods we have em- 
ployed are consistent across the redshift ranges of 
the two samples. 



6.2. Ratios of PAH emission band fluxes 

It has been shown, both experimentally and 
theoretically, that the ratios between the fluxes 
of the 6.2,7.7, 8.6 and 11.3^m PAH emission fea- 
tures are sensitive to the radiation environments 
of the photodissociation regions in which they are 
produced, as well as the sizes of the aromatic 
molecules and the degre e to which they have been 
repro cessed in shocks ( Diamond-Stanic &: Riekd 
2010l and references therein). Since the PAH fea- 
tures are used as a star formation indicator for 
AGN host galaxies, it is important to investigate 
whether the PAH emission band spectra are af- 
fected by the presence of AGN nuclei. 

In an extensive study of the PAH ratios mea- 
sured for a large sample of Seyfert galaxies 
from the Revised Shapley - Ames catalogue (RSA), 
Diamond-Stanic fc Riekd (|2010 h have found that, 
on average, the Seyfert galaxies show 7.7/11.3 
PAH flux ratios that are low compared with local 
star forming galaxies, with some Seyfert galax- 
ies identified as having unusually small ratios 
(7.7/11.3 < 2) ~ lower than measured for any of 
the local s tar forming HII g alaxies in the SINGs 
sample of Smith etal] (|2007i ) . This agrees with 



previous work that also suggested a link between 
low 7. 7/11.3 ratios and the presence of AGN ac- 
tivity (|Smith et al]l2007l : fo'Dowd et al.,,200a) . It 
is therefore interesting to investigate whether the 
radio galaxies in our sample with high EW PAH 
emission — whi ch encompass larger AGN lumi- 
nosities than the iDiamond-Stanic fc Riekl (|2010l ) 
sample — also show low 7.7/11.3 ratios. 

Referring back to Figure [3 the apparent sim- 
ilarity between PAH spectra of the objects with 
high EW PAH detections is striking. Considering 
only the seven objects for which the 6.2, 7.7 and 
11.3/im features can be accurately measured (6 ob- 
jects from the 2Jy sample, and one object from the 
3CRR sample) the mean ratios are: 6.2/11.3 = 
0.7±0.2, 7.7/11.3 = 1.7±0.4, 6.2/7.7 = 0.4±0.1, 
see Table [H It is notable that the 7.7/11.3 ra- 
tios measured for the radio galaxies fall at the 
lower end of the range measured for RSA Seyfert 
galaxies and the Spitzer Infrared Nearby Galaxies 
Survey (SINGs) HII galaxies; two radio galaxies 
(PKS0023-16, 3C293) have 7.7/11.3 ratios that are 
comparable with the most extreme Seyfer t galax - 
ies measured bv .Diamond-Stanic fc Rieke ( 2010[ ). 
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Based on the strong correlation they found be- 
tween 7/7/n^3_j;atios_a£d_^ emis- 



sion, 



Diamond-Stanic fc Riekj (|2010l ) suggested 



that the low 7.7/11.3 may be due to the reprocess- 
ing of the PAH-emitting molecules in shocks, since 
the rotational H2S(3) feature is enhan ced in rela- 
tively hot, shocked regions of the ISM ([Ogle et al 



eiy 

20ia iGuillard et al.ll2009D . Although our sam- 
ple is too small for a thorough statistical investi- 
gation, it is notable in that the objects with the 
smallest 7.7/11.3 ratio - PKS0023-26 and 3C293 - 
also have the highest ratios of i72('S'3)/(11.3-|-7.7) 
amongst all the radio galaxies in the combined 2Jy 
and 3CRR sample (see Table [6|). 

In summary, our results on the PAH ratios for 
the radio galaxies with RSFA are consistent with 
those obtai ned for the most extreme RSA S eyferts 
studied by Diamond-Stanic fc Rieke ( 201Cll ). rein- 
forcing the view that the 6.6, 7.7 and 8.6/zm fea- 
tures are suppressed relative to the 11.3/im feature 
in the environments of AGN, perhaps as a conse- 
quence of reprocessing of the aromatic molecules 
in shocks. Moreover, foll o wing the arguments in 
Diamond-Stanic fc Riekj (|201(i ). we would also 
expect the 11.3/im feature to be a good indicator 
of RSFA in radio galaxies, as it is in other classes 
of AGN. 

6.3. The origin of the far-IR emission 

Based on the correlations between MFIR and 
[OIII] luminosities, we previously concluded that 
the cool, far-IR emitting dust is predominantly 
heated by AGN illumination for the majority 
of radio galaxies that show no independent ev- 
idence for star formation activity (T07, D09, 
DIO). Under the assumption that that 11.3/im 
PAH luminosity is a measure of the overall level 
of star formation activity, if the far-IR con- 
tinuum emission is dominated by AGN-heated 
dust, then we expect the ratio of the PAH lu- 
minosities to the far-IR continuum luminosi- 
ties to be lower than measure d for pure star- 
burst objects. In this vein, ISchweitzer et al 



e gions c 



found that the L-j^-j^rn/vLv^Q^rn ratio 
for PG quasars with detected PAH features 
(i7.7^m/i^it'60Mm=0-0110='=0-002) was similar to 
that of a sample of 12 starburst-dominated 
ULIRGS (L7.7»,m/t^ivi^finum=0.0130 ±0.002), a re- 
sult t hat was also discussed by IVeilleux et al. 
()2009l) . The authors conclude that the similar 
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Fig. 14. — : Plot of Lii.3^m/i^ii/60/im vs 
vL^dOpLTa.. The 60/im fluxes for the starburst 
galaxies were obtained from the IRAS catalogue, 
and the 2Jy and 3CRR sample 70/im fluxes were 
interpolated to 60/im using the 70/24/im spectral 
index from the Spitzer data. The IRS data for the 
starburst galaxies and ULIRGs were downloaded 
from the Spitzer Heritage Archive and re-reduced 
and analyzed in an identical way as described for 
the radio galaxies in this paper. In addition, to 
compensate for flux loss in the IRS spectra due 
to the extended structure of the starburst galaxies 
outside of the slit, a correcti on factor was applied 
as discussed in lBrandl etahl ([2006). The ULIRGs 
sample comprises 15 objects from the sample pre- 
sented in Veilleux with an average AGN contribu- 
tion of <20%. The archetypal ULIRG ARP220 is 
marked with a red diamond. 



ratios argue in favor of a starburst origin for the 
far-IR emission of PG quasars. In cont rast, based 
on similar techniques, IShi et al.l (|2007l ) concluded 
that the contribution of starburst heating to the 
far-IR continuum is relatively minor (~24%) for 
most radio-loud AGN. 

The problem with using the LpAn/Lfar-iR ra- 
tio as a direct indicator of the degree of starburst 
heating of the far-IR emitting dust is that the 
scatter in this ratio is large, even for starburst 
dominated objects. Results from SINGs have in- 
dicated that the use of PAH emission for esti- 
mating star formatio n rates may be uncertain by 



20071) . 



20051 : ICalzetti et all 



a factor of ~10 (Da le et al. 

This conclusion is also supported by 
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Table 6 
PAH Ratios 



Name 


6.2/11.3 


7.7/11.3 


6.2/7.7 


H2(S3)/(11.3 + 7.7) 


0023-26 


0.8 ± 0.1 


1.2 ± 0.1 


0.7 ± 0.1 


0.113 ± 0.01 


0347+05 


0.7 ± 0.1 


1.8 ± 0.1 


0.4 ± 0.1 


0.008 ± 0.01 


0915-11 


0.8 ± 0.1 


1.8 ± 0.2 


0.4 ± 0.1 


0.005 ± 0.01 


1733-56 


1.0 ± 0.1 


1.9 ± 0.2 


0.5 ± 0.1 


0.008 ± 0.03 


2135-20 


0.5 ± 0.1 


1.9 ± 0.2 


0.3 ± 0.1 


0.018 ± 0.02 


2314+03 


0.8 ± 0.1 


2.4 ± 0.2 


0.3 ± 0.1 


0.011 ± 0.01 


3C293 


0.4 ± 0.1 


1.2 ± 0.1 


0.3 ± 0.1 


0.156 ± 0.01 



Note. — Table presenting the ratios of PAH and H2 fluxes for objects with high 
EW detections of PAH. Sec Section [6?2] of the discussion. This table does not 
include 3C285 and 3C305 that have upper limits for their 6.2 and 7.7/im PAH 
fluxes. 



Tacconi-Garman et al. (I2005I) who find strong 
variations in the PAH to stellar continuum ration 



in nea rby starburst galaxies, as well as lHaas et al 



( 20091 ) who find a range in the L/ir/.Z^6.3pm of a 
factor of 10 for a sample of PAH selected galaxies 
from the ISOCAM parallel survey. In Figure [T3] 
we plot _Lii.3^,„/i/L^60/im for the 2Jy and 3CRR 
objects from the current study, along with com- 
parison samples of LIRGs and starb urst objects of 
simila r far-IR luminosity taken from ISrandl et al.l 
( 20061 ). and star burst-dominated ULI RGs taken 
from the study of IVeilleux et al. I (I2OO9I) . 

Consistent with the above discussion it is no- 
table that the scatter in the Lii.3^m/i/Lj,60/im 
ratio is large, both for the radio galaxies we 
have identified as RSFA objects in our sam- 
ple, a nd for the starburs t galaxies in the sam- 
ple of iBrandl et al. In addition, and pei- 
haps counter-intuitively, the starburst-dominated 
ULIRGs tend to the lowest values of Lii.3^,„/i/Lj,60^m. 
An explanation for this difference might be that a 
significant fraction of the starburst activity in the 
ULIRGs is so deeply embedded and extinguished 
by dust that it is not detected using the short 
wavelength PAH features (see Veilleux et al. 2009 
for further discussion). 

Given the high degree of scatter and large num- 
ber of relatively high upper limits in the plot, it 
is difficult to make strong conclusions about the 
heating mechanism for the far-IR emitting dust 
for most of the radio galaxies in our sample based 
on this evidence alone. However, the upper limits 
for some radio galaxies fall well below the majority 
of values measured for the Brandl et al. LIRG and 



starburst galaxies {Lii,^^^/ ^L^&Qy,ui < 5x lO^^^), 
suggesting that AGN illumination makes a large 
contribution to the heating of the cool dust in 
these objects, supporting our previous conclusions 
(T07, D09). It is noteworthy that many of the 
objects with LpAHii.Stim/i^L^ 60^m < 5 x 10"'^ 
comprise the low-z NLRG with high S/N spec- 
tra (e.g. PKS1559+02, PKS1949-h02, PKS2250- 
41, PKS2356-61) as weU as low-z BLRG (3C382, 
3C390.3) with the very highest S/N spectra. This 
suggests that many of the objects with low S/N 
spectra (mainly high-z objects) might fall in this 
region of the diagram if they had better S /N data. 

These results are relevant to the interpretation 
of the far-IR emission from AGN in the distant 
Universe. High redshift galaxie s that show evi- 



dence for AGN a ctivity at X-ray (lAlexander et al 



2OO 5I). optical (IPriddev et all I2OO3I) and radio 



(Archibald et al.ll200U IWillott et al.l 12002 ) wave 



lengths have been detected with deep surveys at 
sub-mm wavelengths which are sensitive to the 
redshifted far-IR emission from cool dust. Be- 
cause investigations of starbursts galaxies in the 
local universe have revealed prodigious far-IR ra- 
diation, the sub-mm emission has been interpreted 
as being due to cold dust heated by strong star 
formation activity. Our previous work on the 2Jy 
sample (T07, D09, DIO) has shown that this inter- 
pretation is not unique because the dust emitting 
the thermal mid- to far-infrared (MFIR) radiation 
can be heated by either starbursts or by AGN. Fig- 
ure [14] provides evidence that, in many powerful 
radio galaxies, far-infrared emission is dominated 
by AGN and not starburst heating. Therefore, our 



26 



results suggest that the interpretation of sub-mm 
emission from distant radio galaxies as originating 
purely from star formation should be treated with 
caution. 

6.4. Compact radio sources 

Compact steep spectrum (CSS) and gigahertz 
peaked sources (GPS) are characterized by the 
small sizes of their radio structures (D < 30 
kpc and <lkpc respectively) and steep high fre- 
quency radio spectra. Mounting evidence sup- 
ports the idea that CSS/GPS objects represent a 
young phase of more extended powerful radio- 
loud AGN (|Fanti et al.l Il995l iReadhead et al.l 
19961 ). In particular, the measured proper mo- 
tions of the radio components in some CSS objects 
that imply very young ag es (<3xl0'^ years; e.g. 



Polatidis fc Conwavll2003D . 

The combined 2Jy and 3CRR sample includes 
7 CSS sources and 1 GPS source (PKS1934-63). 
Three of these objects have high equivalent with 
PAH detections (PKS0023-26, 3C305, PKS2135- 
20) and three additional objects have been iden- 
tified with low EW PAH detections (PKS1151-34, 
PKS1814-63, PKS1934-63). The rate of CSS/GPS 
objects in the 2Jy-f 3CRR sample with PAH detec- 
tions is therefore 75% (6/8) — much higher than 
the detection rate for the extended radio sources 
in the 2Jy-|-3CRR sample with good IRS spectra 
(21%: 10/48). Int erestingly, similar re sults have 
been obtained bv IWillett et al.l (|2010f ). who de- 
tect PAH emission in 7/8 (88%) of the objects 
in their heterogeneous sample of nearby compact 
symmetric object^ (z < 0.26). This evidence for 
an enhanced degree of star formation in the host 
galaxies of compact radio sources is also consis- 
tent with the fact that the compact sources are 
disproportionately represented in the (small) sub- 
set of radio galaxies that show strong evidence for 
young stellar populatio ns in their optical spectra 



( Tadhunter et al.ll201lD . 



At first sight it might be thought that the rel- 
atively high incidence of starburst activity among 
the compact radio sources is connected with the 
relative youthfulness of such sources: if the ra- 



^The CSOs represent a subset of the compact radio source 
population, with relatively symmetric radio structures; 
most of the compact (CSS, GPS) sources in the 2Jy-|-3CRR 
sample would also be classified as CSOs. 



dio source activity is triggered close to the peak 
of the starburst associated with the triggering 
events, then the compact sources might be ex- 
pected to show enhanced star formation activ- 
ity compared to the more extended radio sources 
that are observed at later epochs. However, the 
problem with this explana t ion is that both models 
(jMihos fc HernquistI 11995 toi Ma tteo et al.ll2007l) 



and direct observations fiRodrig uez Zaurm et al 
20101) suggest that the main starburst episodes 



associated with major galaxy mergers are rela- 
tively long-lived: ~10 — 100 Myr - much longer 
than the expected time period over which a radio 
source would appear as a compact radio source, 
but similar to expected lifetimes of extended ra- 
dio sources. Moreover, the spectroscopic signa- 
tures of the merger-induced starburst would be 
detectable in optical spectra for several hundred 
Myr after the starburst. Therefore we would not 
expect to find a major difference between the rate 
of detection of starburst signatures in compact and 
extended radio sources, given the similarities be- 
tween the expected timescales of (extended) radio 
sources and starburst activity. 



Alternativelv. lTadhunter et al.l(|2011l ) have sug- 
gested that the association between starburst ac- 
tivity and compact radio sources may be a conse- 
quence of an observational selection effect: young 
radio sources triggered in major gas-rich mergers, 
that also lead to starbursts, will have their ra- 
dio emission enhanced by jet/cloud interactions as 
their radio jets plough through the dense merger 
debris. This will lead to such sources being pref- 
erentially selected in flux limited radio samples 
such as the 2Jy and 3CRR, boosting their num- 
bers relative to the radio sources that are not trig- 
gered in major gas-rich mergers; the compact radio 
sources triggered in this way will have a higher ra- 
tio of radio power to intrinsic jet and AGN power 
than more extended radio sources triggered in less 
dense environments. As these compact sources 
expand further into the more tenuous haloes of 
the host galaxies their radio fluxes will decline 
and they will drop out of the flux limited sam- 
ples. We note that this selection effect is consis- 
tent with the fact that the compact radio sources 
tend to fall below the correlations between [OIII] 
and 24/Ltm emission and radio power, as expected 
if the radio flux is boosted relative to the intrin- 



sic AGN power (iMorganti et al.lll997t iHolt et al 
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20091 : Mo rganti et al. 2011). It m ay also help ex- 
plain why lO 'Pea fc BaumI (|l997t ) found that com- 
pact radio sources are more common in flux lim- 
ited radio samples than would be expected on the 
basis of a simple extrapolation of the radio size vs. 
number relation for more extended radio sourcej^ 

6.5. Consequences for the triggering and 
evolution of radio source host galaxies 

In parallel work, we have found that at least 
85% of the objects in the 2Jy sample show 
morphological evidence f or recent galaxy in- 



teracti ons and mergers ([Ramos Almeida et al 
2011bl la). However, the morphologies and stel- 



lar population properties displayed by the ra- 
dio source host galaxies are diverse and do not 
all correspond to a single phase of a particular 
type o f galaxy interaction ( Ramos Almeida et al 



2011bllal : iTadhunter et alll201ll) . While some ra- 
dio sources have indeed clearly been triggered at 
the peaks of major gas-rich mergers, others are 
seen at a relatively late post-merger phases well 
after the coalescence of the merging nuclei; and 
a significant subset show tidal links with well- 
separated companion galaxies, suggesting that 
they have been triggered after the first pass of 
the interacting galaxies, but well before the final 
merger (if any). These results indicate that it 
is possible for the gas inflows associated with a 
range of stages and types of galaxy interactions 
to be sufficient to trigger powerful radio galax- 
ies, many of which are associated with AGN of 
quasar-like luminosity {Lbol > 10'^® W). Over- 
all, these results are inconsistent with many re- 
cent hydrodynamical simulations which predict 
that the most powerful, quasar-like AGN activity 
is only triggered around the time of coalescence of 
the merging nuclei in major, gas-rich mergers (e.g. 
di Matteo et al.ll2005l ). or with a short time delay 
(|Hopkinsll201l[) . Our mid-IR results on the PAH 
features strongly reinforce these conclusions: in 
most objects we do not detect the levels of RSFA 



l^Whereas the extended radio sources (D > 10 kpc) show a 
rapid increase in number within bins of increasing Hnear 
size for the radio sources, the compact radio sources {D < 
10 kpc) show a constant number with size — there are 
far more compact sources than expected on the basis of a 
simple extrapolation of t he number vs size rela tionship for 
the extended sources fsee lO'Dea fc Baum|[l997l. Figure 10) 
towards small sizes. 



that would be expected if all powerful radio galax- 
ies were triggered at the peaks of major, gas-rich 
merger|3 We suggest that the current hydrody- 
namical simulations do not adequately capture 
the physics of the gas flows across the variety of 
interaction types and/or stages that lead to the 
triggering of radio-loud AGN. 

It is also important to consider whether the neg- 
ative feedback effect of the AGN and jet activity 
might be responsible for the low levels of RSFA 
detected in radio galaxies i.e. whether the AGN 
and jets in these powerful radio galaxies input suf- 
ficient energy into the ISM of the host galaxies 
to suppress some of the star formation associated 
with the triggering events. Because star formation 
associated with advanced mergers is often concen- 
trated near the nucleus, the energy input of the 
AGN may have a major impact on these star for- 
mation regions. Given that we have found evi- 
dence for RSFA in up to 35% of our samples, it 
appears that the AGN cannot entirely shut down 
star formation in powerful radio galaxies. On the 
other, we note that two of the objects in this inves- 
tigation (3C293, 3C326) have been associated with 
AGN feedback processes through detailed observa- 
tions of the mo lecular gas content. The study of 
iNesvadba et al. found evidence that these 

objects have star formation rates lower than ex- 
pected for the amount of dense gas available for 
star formation. 

The feedback effect of the AGN may be sig- 
nificant in some objects, however, the lack of ev- 
idence for RSFA in the 2Jy and 3CRR samples 
may simply be due to the fact that we are ob- 
serving the sources a long time before or after any 
main merger-induced starburst. Alternatively, the 
triggering galaxy interactions may not be of the 
types that are associated with major star forma- 
tion episodes (e.g. they may be relatively minor 
or gas-poor). Clearly, in order to develop our un- 
derstanding of the triggering events further, we re- 
quire information about the total reservoir of cool 
ISM in the host galaxies, as well as the details of 
the circum-nuclear gas kinematics. 



^^Note that we find no clear link between those objects 
we have identified in this study with evidence for RSFA 
and any one particular stage or type of interaction (i.e. 
with signatures of pre-coalescence, coalescence and post- 
coalescence) that have been identified in the host galaxy 
morphologies. 



28 



7. Conclusions 

We have presented Spitzer/IRS spectroscopy 
for complete samples of 2 Jy and 3CRR radio-loud 
AGN, detecting 93% and 79% of the objects in the 
two samples respectively. Our analysis of the spec- 
tra reveals strong RSFA-tracing PAH features in 
only a minority of the objects from the two sam- 
ples (16%) that have good IRS spectra. Com- 
bining this result with optical continuum spec- 
troscopy, mid- to far-IR (MFIR) color and far-IR 
excess diagnostics, wc find that only 35% of ob- 
jects in the combined 2Jy and 3CRR sample show 
any evidence for RSFA at optical and/or MFIR 
wavelengths. This result argues strongly against 
the idea that there is a close link between starbursts 
and powerful radio-loud AGN, reinforcing the view 
that only a minority are triggered at the peaks of 
star formation activity in major, gas-rich mergers. 
The PAH emission also allows us to test whether 
there exists a substantial proportion of RSFA that 
is obscured by dust that was undetected at opti- 
cal wavelengths. Although we do find significant 
PAH emission in 4 objects that were not previ- 
ously identified as hosting RSFA on the basis of 
other star formation indicators, this does not sub- 
stantially change the statistics for the overall rate 
of detection of RSFA in radio galaxies from our 
previous investigations. In addition, we find that 
compact radio sources show a significantly higher 
incidence of RSFA, which cannot be readily ex- 
plained by the youthful nature of these objects. 
We suggest that, for radio selected samples, there 
may be a bias towards the selection of compact 
radio sources that are triggered in gas-rich envi- 
ronments. 
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Table 7 
Ancillary data 



Name z Sii.3(W/m^) S6o(W/Hz) 



ULIRGs 



F12112+0305 





0733 


8 


5E- 


16 


2 


5E+26 


F17207-0014 





0428 


3 


5E- 


15 


2 


9E+26 


F20414-1651 





0871 


3 


6E- 


16 


1 


8E+26 


F13335-2612 





1250 


3 


5E- 


16 


1 


3E+26 


F03250+1606 





1290 


3 


5E- 


16 


1 


4E+26 


F10565+2448 





0431 


2 


8E- 


15 


1 


lE+26 


F23234+0946 





1279 


1 


7E- 


16 


1 


7E+26 


F14348-1447 





0830 


7 


6E- 


16 


2 


6E+26 


F22491-1808 





0778 


3 


lE- 


16 


1 


7E+26 


F16333+4630 





1910 


3 


OE- 


16 


3 


4E+26 


Arp 220 





0181 


3 


5E- 


15 


1 


6E+26 


F14197+0813 





1310 


1 


lE- 


16 


1 


lE+26 


F05024-1941 





1920 


1 


lE- 


16 


2 


6E+26 


F01494-1845 





1580 


3 


5E- 


16 


2 


OE+26 


F00482-2721 





1292 


9 


3E- 


17 


1 


lE+26 



LIRGS and staibuists 



IC 342... 





0001 


1 


8E- 


14 


8 


3E+21 


Mrk 52... 





0071 


2 


lE- 


15 


1 


lE+24 


Mrk 266... 





0155 


2 


7E- 


15 


7 


9E+24 


NGC 520... 





0076 


2 


9E- 


14 


8 


lE+24 


NGC 660... 





0028 


4 


9E- 


14 


2 


3E+24 


NGC 1097... 





0042 


4 


3E- 


14 


4 


2E+24 


NGC 1222... 





0082 


6 


5E- 


15 


3 


8E+24 


NGC 1365... 





0055 


3 


OE- 


14 


1 


2E+25 


NGC 1614... 





0159 


1 


6E- 


14 


3 


6E+25 


NGC 2146... 





0030 


1 


6E- 


13 


5 


7E+24 


NGC 2623... 





0185 


4 


2E- 


15 


3 


7E+25 


NGC 3256... 





0094 


4 


lE- 


14 


4 


OE+25 


NGC 3310... 





0033 


1 


6E- 


14 


1 


7E+24 


NGC 3556... 





0023 


3 


OE- 


14 


7 


7E+23 


NGC 3628... 





0028 


3 


6E- 


14 


1 


9E+24 


NGC 4088... 





0025 


1 


9E- 


14 


7 


5E+23 


NGC 4194... 





0083 


1 


6E- 


14 


7 


lE+24 


NGC 4676... 





0220 


3 


7E- 


15 


5 


9E+24 


NGC 4818... 





0036 


8 


9E- 


15 


1 


lE+24 


NGC 4945... 





0019 


1 


9E- 


13 


9 


7E+24 


NGC 7252... 





0160 


6 


OE- 


15 


4 


6E+24 


NGC 7714... 





0093 


7 


OE- 


15 


4 


3E+24 



Note. — Tabic presenting the ancillary data from Figure 
1141 for ULIRG, LIRG and starburst objects. 
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